Phase diagrams of Fermi gases in a trap with mass and population imbalances at finite temperature 
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The pairing and superfluid phenomena in a two-component Fermi gas can be strongly affected by the pop- 
ulation and mass imbalances. Here we present phase diagrams of atomic Fermi gases as they undergo BCS- 
Bose-Einstein condensation (BEC) crossover with population and mass imbalances, using a pairing fluctuation 
theory. We focus on the finite temperature and trap effects, with an emphasis on the mixture of ®Li and '"^K 
atoms. We show that there exist exotic types of phase separation in the BEC regime as well as sandwich-like 
shell structures at low temperature with superfluid or pseudogapped normal state in the central shell in the BCS 
and unitary regimes, especially when the light species is the majority. Such a sandwich-like shell structure ap- 
pear when the mass imbalance increases beyond certain threshold. Our result is relevant to future experiments 
on the ''Li~'*"K mixture and possibly other Fermi-Fermi mixtures. 

PACS numbers: 03.75. Ss,03.75.Hh,67.85.Pq,74.25.Dw 



Ultracold Fermi gases provide an excellent model system 
for studying condensed matter physics, owing to the vari- 
ous experimentally tunable parameters. Using a Feshbach 
resonance, a two-component Fermi gas of equal spin mix- 
ture of atoms of equal mass exhibits a perfect crossover 
from Bardeen-Cooper-Schrieffer (BCS) type of superfluid- 
ity to Bose-Einstein condensation (BEC) jr-^, which has 
been studied intensively both in experiment |4-8] and theory 
I9I, Toll. There have also been a great deal of experimental 
11, 12] and theoretical ][l3-18] studies on equal-mass sys- 



tems with population imbalance. Population imbalance adds 
a new dimension to the phase diagrams, leading to phase 
separation 119], Sarma superfluid 120] and possibly Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) states 112 ill . Mass imbal- 
ance, i.e., pairing of different mass atoms, will further enrich 
the physics. Indeed, there have been world wide efforts on 
the study of Fermi-Fermi mixture of different species. Over 
the past several years, Feshbach resonances between differ- 
ent species of fermionic atoms, e.g., ^Li and ^''K, have been 
found and studied 1.22l424ll . although it remains to achieve su- 
perfluidity experimentally. There have been some theoretical 
studies in this aspect, e.g., on the strong attraction limit at 
zero temperature T li25il . few -body physics [26] or the po- 
laron physics i27 [ | 28ll as well as thermodynamics of a high T 
normal mixture [29] . There have also been studies on phase 
diagrams, which, however, are mostly restricted to zero tem- 
perature using a mean field theory either in a homogeneous 



Fermi gas fsH.'sT] or in a trap |l32-35]. Recently, Guo et al. 

studied the mass imbalanced Fermi gases at finite temper- 
atures. Due to technical complexity, this study was restricted 
to homogeneous systems only. In order to address various ex- 
periments, which are always done at finite T and in a trap, it is 
important to take into account the trap and finite temperature 
effects simultaneously. 

In this paper, we consider a two-species Fermi-Fermi mix- 
ture with a short-range s-wave pairing interaction in a 3D 



isotropic harmonic trap at finite temperature. We emphasize 
on the interplay between the finite T ij^ 37] and trap effects 
while the mass ratio m^/mj, and population imbalance (or 
"spin polarization") p = (iV-|- — Ni)/{N^ + N{) (as well as 
the ratio / between the trapping frequencies) are varied 
within a pairing fluctuation theory, where spin index a =tj i 
refers to the heavy and light species, respectively. In order 
to address experiments, we pay special attention to the ^Li- 
mixture, while keeping in mind other possible mixtures 
such as '^Li-i'^^Yb and i^^Yb-^'^^Yb. We will present our 
theoretical findings in terms of representative phase diagrams 
in the T-p plane, as long with typical gap and density pro- 
files, throughout the BCS-BEC crossover One special fea- 
ture here is the emergence of wide-spread pseudogap phenom- 
ena at finite T. For high enough mass imbalance, our result 
shows that, in the T-p phase diagrams, three-shell sandwich- 
like spatial structures occupy a large region both at unitarity 
and in the BCS regime, including sandwiched phase sepa- 
ration, sandwiched Sarma, sandwiched polarized pseudogap 
states with increasing T. In the BEC regime, there are exotic 
"inverted" phase separations with a normal Fermi gas core in 
the trap center surrounded by paired (superfluid or pseudogap) 
state in the outer shell. Our results provide excellent predic- 
tions for future experiments. 

Except for slightly different notations, our formalism is a 
combination of that used in Refs. lfl6l, [36ll . where a single- 
channel Hamiltonian is used to describe the Fermi gas, with 
a local density approximation (LDA) for addressing the trap 
inhomogeneity. The heavy (cr and light {a =\) species 
have dispersion ^y^r = eko- — l^^a = h!^ /2m„ — fic, and bare 
fermion Green's function Gq^{K) = za;„ — ^ko-, where /io- 
is the chemical potential, w„ the fermionic Matsubara fre- 
quency. As in Ref. ^3^, we take h = ks = I and use the 
four vector notation, e.g., K = (w„,k), J2k = ^Z^nSk' 
etc. At finite T, the self-energy Yi^iK) contains two parts, 
Eg. (if) = Ssc,(t(^) + "^yg.aij^), where the condensate 
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contribution T,sc,a{K) = —Al^Goai—K) vanishes above 
Tc, and the finite momentum pair contribution Epg ct(^) = 
^Qtpg{Q)Goa-{Q — K) persists down to T = 0. Here 
a = —a, and the T-matrix tpg(ff) = g/[l + gxiQ)] repre- 
sents an infinite series of particle-particle scattering processes, 
with a short-range interaction strength .g < and the pair sus- 
ceptibility xiQ) = Y.K.aGoa{Q - 'k)Gs[K)/2. At and 
below Tc, we have Y.pg,„{K) « Y.Qhg{Q)GQs{-K) = 
— Ap^Goct(— ^) + <^S, which defines a pseudogap Ap^ via 
= — TliQ tpgiQ)- Ignoring the less important incoherent 
term (JS, we obtain ^^(X) = —A'^Goai—K) in the sim^ 
pie BCS form, where A^ = A^ 
Green's function is given by 
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A^g. Therefore, the full 



GaiK) = 
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where 

V1TA2, and Et,„ ^Ek + Ck., ^k = (^t + Cu)/^^ Ck. - 
(^k<T - Ck<T)/2. With rv = Y.K Ga{K), n ^ + 71^ and 
6n = n-^- — n^, the number equations read 
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where the average Fermi function f{x) = /(^ + Ckcr)/2. 

At T < Tc The Thouless criterion leads to the gap equation 
9^^ + X(0) = 0. For T > it is amended by g-^ + x(0) = 
ZjjLp, where the effective pair chemical potential /ip and the 
coefficient Z can be determined from the Taylor expansion 
of the inverse T-matrix [9], tpg{Q) = Z{iQi - Qq), with 
rtq = /2M* — jip being the pair dispersion, and M* the 
effective pair mass. Thus the gap equation reads 
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1 - 2/(Sk) 



L2ek 



(4) 



with /ip at T < Tc. Here g is regularized by = 
mr/2Tia — l/2ek' where a is the s-wave scattering length, 
rrir = m^mi/ {m^ + mi) the reduced mass, and = 
k'^/Anir = Ck + with /i = (^^ + ^;)/2. 

The T-matrix expansion leads to the pseudogap equation 



(5) 



where b{x) is the Bose distribution function. As in Ref. |l3f 
we impose a cutoff qc on the summation such that pairs with 
q > Qc may decay into the particle-particle continuum. 

In a harmonic trap, the LDA approximation imposes that 
the local /Ict (r) — /io-lO) — ^TOo-w^^^, with trap frequency ojo-. 
We have the total particle number N ~ J d'^r n{r) and the 
number difference 6N = N^—N^ = pN = J d^r Sn{r). The 



Fermi energy = {SNy/^w^ = kj,/2r 
Tp is defined as that for an unpolarized, noninteracting Fermi 
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Figure 1. (Color online) T~p Phase diagram of ^Li-'^'K mixture m 
a harmonic trap at unitarity, with u)-^ = ^i- The solid lines separate 
different phases and the (red) dashed line between the (Sandwiched) 
PG and the normal phases is approximated by mean field calcula- 
tions. Here "PG" and "PS" indicate pseudogapped normal state and 
phase separation, respectively. 



gas with the same total number N and trap frequency w^, 
with TO = (to-j- + mi)/2, the average mass. Here Rtf is 
the Thomas-Fermi radius, and the species dependent R^p — 

^2{6Nay/^/maOJa. 

For a homogeneous system, Eqs. (|2]l-(|5]l form a closed set 
and can be used to solve for Tc, as well as /i. A, A^c and Apg 
for given T. However, in a harmonic trap, this set of equations 
have to be solved at each given radius and then subject to the 
total trap-integrated particle number constraints. 

In order to find the stable states in a trap, we compare the 
(local) thermodynamical potential ils (per unit volume) in a 
superfluid or pseudogap state with its normal Fermi gas coun- 
terpart, = -TJ2k.^ln.{l + e^^^"/'^). The paired state 
rts consists of contributions from fermionic excitations flp 
and noncondensed pairs ils- 



nF = -—+ E(ek - E,) tE ln(l + e-^-/^), 



(6) 



fiB = Z/ipA2^ + TEln(l 



The stable states should have a lower value. When fts > 
at certain radii r, phase separation takes place. Note that the 
Hb term is absent in simple mean-field calculations. 

Shown in Fig.[T]is the calculated T-p phase diagram at uni- 
tarity 1/kpa = with = uj^ and mass ratio mf/mi = 
40:6, as is appropriate for the ^Li-^"K mixture. Correspond- 
ing representative density and gap profiles are shown in Fig.|2] 
It can be seen that phase separation (PS) and sandwiched PS 
occupy the lowest T part. In particular, when the light species 
dominates the population, we have a regular phase separation 
with an equal population superfluid core in the trap center. 
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Figure 2. (Color online) Representative density (main panels) and 
gap (insets) profiles for p — —0.85 and 0.25 at unitarity: (a)-(c) 
T/Tf = 0.2, 0.1, 0.02; (d)-(f) T/Tf = 0.15, 0.05, 0.01, corre- 
sponding to each phase in Fig.[T] The temperatures in (c) and (f) are 
essentially zero. The (black) solid and (red) dashed lines in the main 
panels are for ®Li (nj_, light species) and ^''K (n^, heavy species), re- 
spectively. Plotted in the insets are the total gap A (black solid lines) 
and order parameter Asc (blue dashed lines). Here the insets share 
the same x-axes, and the densities Ua- are in units of Up = k% /Stt^. 



surrounded by the majority light atoms in the outer shell, sim- 
ilar to that seen in the equal-mass case lfl6ll . Here the gap A, 
the order parameter Age, and the density of minority atoms 
ni jump to zero across the interface, as shown in Fig. |2c). 
However, except for the light atom dominated case, a three 
shell structure appears, which we refer to as sandwiched PS. 
Typical density and gap profiles are shown in Fig.|2|f) (for es- 
sentially zero T), where an equal population superfluid exists 
only at intermediate radii, whereas the light and heavy atoms 
dominate the outer and inner shells, respectively. In fact, at 
zero T, the heavy atoms are absent in the outer shell. The 
gaps and density profiles exhibit first order jumps at both in- 
terfaces. Our findings at the lowest T are consistent with the 
earlier works at zero temperature 11321 13511 . The primary rea- 



son for this three shell structure to occur is that R^p <IC R^p 
so that at the trap center It is known that at low 

T, population imbalance tends to break pairing OlSll . There- 
fore, only at certain intermediate radii where Uf « nj, can a 
superfluid exist. As p decreases from 1 towards -1, these radii 
move from the trap edge near Rj,p towards the trap center, 
until the inner shell of normal mixture shrinks to zero around 
p = -0.71. 



Previous study OlSIl demonstrates that in a homogeneous 
Fermi gas of equal mass at unitarity, a population imbalanced 
Sarma superfluid ||20|] is stable only at intermediate temper- 
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Figure 3. (Color online) Phase diagram of ''Li-'"'K as in Fig.[T]but 
for (a) l/kpa = -0.5 and (b) -^0.5. Here "PS-SF" and "PS-PG" 
represent phase separated superfluid and phase separated pseudogap 
phases, respectively, with a normal gas core surrounded by an outer 
shell of superfluid or pseudogapped normal mixture. 



atures. Here we find this remains true for an unequal mass 
mixture. In Fig. [H the PS phase becomes a Sarma phase at 
intermediate T, where population imbalance penetrates into 
the inner superfluid core so that the first order jumps of the 
gap, the order parameter and the minority density at the in- 
terface disappear This is similar to the Sarma state found 
in Ref. d. Correspondingly, the sandwiched PS phase be- 
comes sandwiched Sarma phase. As seen from the p = —0.85 
and 0.25 cases shown in Figs.|2|b) and|2|e), respectively. A, 
Asc and vanish continuously at the interface (of the larger 
radius). The difference between A and Asc defines the pres- 
ence of the pseudogap (not shown) Apg — ^A^ — Age- 

As the temperature increases further, the superfluid region 
disappears so that the Sarma phase becomes a polarized pseu- 
dogap (PG) phase, where a finite pseudogap exist in the in- 
ner core without superfluidity. In a similar fashion, the sand- 
wiched Sarma phase evolves into a sandwiched PG phase. 
Representative density and gap profiles are shown in Figs.|2a) 
and Ed) for (p, T/Tp) = (-0.85, 0.2) and (0.25, 0.15), re- 
spectively. Finally, at very high T we have a normal phase. 
Note that the separation between the normal and the PG or 
sandwiched PG phases is a crossover rather than a phase tran- 
sition, and thus is determined approximately using the BCS 
mean-field theory. 

Next, we present in Fig.[3]the phase diagrams at (a) \lkFa=- 
0.5 and (b) 0.5, similar to Fig. [T] but in the (near-)BCS and 
(near-)BEC regimes, respectively. For the BCS case, except 
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Figure 4. (Color online) Typical density and gap profiles for BCS 
(left column, l/kpa — —0.5) and BEC (right column, l/kpa = 
0.5) regimes. The convention is the same as in Fig.|2] Panels (a) and 
(b) are for p = —0.25 at T = O.OITf ~ 0, corresponding to the 
sandwiched PS phase in Fig.[3ja). Panels (c) and (d) plot the density 
and gap (insets) distributions forp = 0.25 at T = 0.2Tf and O.ITf, 
representing the PS-PG and PS-SF phases in Fig.|3tb), respectively. 



for the high T normal phase, the phase diagram is essentially 
occupied by three-shell sandwich-like structures. The cen- 
tral shell is an unpolarized BCS superfluid at the lowest T, 
spin polarized Sarma superfluid at intermediate T, and a po- 
larized pseudogapped normal state at slightly higher T. In the 
sandwiched PS phase, the outer shell is a normal mixture with 
light species (^Li) in excess, surrounded by a single compo- 
nent normal gas of the light atoms at the trap edge. This is dif- 
ferent from the sandwiched PS at unitarity [Fig.|2jf)], where 
the outer shell contains the light atoms only at T = 0. Typical 
gap and density profiles (at p = —0.25 near T = 0) are shown 
in Figs. |4ja) and|4jb), respectively. In comparison with the 
unitary case, one can see that as the pairing strength decreases 
from unitarity, the sandwiched PS phase in the phase diagram 
expands and gradually squeezes out the PS phase completely. 
The evolution of the various phases with increasing tempera- 
ture is similar to their unitary counterparts, except that now the 
sandwiched PG phase occupies a very slim region, reflecting 
a much weaker pseudogap effect in the BCS regime. 

The phase diagram for the BEC case in Fig. |3lb) is rather 
different. First, for p < 0, where the light species is the ma- 
jority, a Sarma superfluid phase occurs at low T, with essen- 
tially equal population at the trap center. Indeed, polarized 
superfluid becomes stable and phase separation is no longer 
the ground state in the BEC regime il5ill9ll . As T increases, 
the order parameter decreases to zero and the system evolves 
into a polarized pseudogapped normal state. The large area of 
the "pseudogap" phase indicates greatly enhanced pseudogap 
effects in the BEC regime. On the other hand, for (roughly) 
p > 0, where the heavy species dominates, we have an "in- 
verted" phase separated superfluid state at low T, labeled as 
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Figure 5. (Color online) Spatial distributions of various phases in the 
trap at unitarity as a function of mass ratio m^/m^ for p = —0.23 
with (a) T = O.ITf, ui^ = oJi and (b) T = 0.15Tf and uji = 2u^. 
The (green) dotted line indicates the mass ratio for ^Li-^^K. 



"PS-SF", where a normal gas core of the heavy species is sur- 
rounded by a shell of unpolarized superfluid. This should be 
contrasted with the PS phase in the unitary case [Fig. |2|c)], 
where the normal Fermi gas is outside the superfluid core. As 
T increases, a phase separated pseudogap state (labeled "PS- 
PG") appears, where pseudogap exists in the polarized outer 
shell but without superfluidity. This is an exotic new phase, 
which has never been seen or predicted before. Typical den- 
sity and gap (insets) profiles for the PS-PG and PS-SF phases 
are shown in Figs.|4fc) and|4fd), respectively. Possible causes 
for the "inversion" of the phase separation include: (i) For 
p > 0, R}pp of the heavy species becomes close to R^p of 
the light atoms; (2) As the local kpir) decreases with r, the 
outer region is deeper in the BEC regime than the trap center, 
making pairing easier and energetically more favorable at the 
trap edge. When compared with the three-shell structure at 
unitarity, one concludes that as the pairing strength increases 
from unitarity, the outer shell of normal light atoms retreats 
and finally disappears. 

We now turn to the case of variable mass ratio rrif/mi, 
with different cj^/cj^. Plotted in Fig. |5] are spatial distribu- 
tions of possible phases at unitarity as a function of the mass 
ratio atp = -0.23 for (a) (T/Tp, WfM) = (0-1,1) and 
(b) (0.15, 1/2), respectively. For both cases, a sandwich-like 
structure appears as the mass ratio increases beyond 3.7 and 
7.0, respectively. This can be easily understood by looking at 
the corresponding non-interacting density distributions 13211 . 
Starting from m^/mi — 1 and lo^/uji = 1, the majority 
species always has a larger spatial extension so that pairing 



5 



is easier at the trap center. However, as m^/rai becomes 
sufficiently large, i?^^ of the heavy species may become 
smaller than R^p so that the non-interacting density distribu- 
tion curves cross each other at an intermediate radius, where 
pairing is more energetically favorable than elsewhere so that 
a three-shell structure appears at low T, as shown in Fig.lSja). 
Note that BCS pairing requires a match of (mass independent) 
k^p = (67r^ncr)^/'' (locally) between the two species. Thus 
the mass ratio changes the position of the crossing point via 
changing R^p- 

Tuning w^/w^ can also change the density crossing point 
since R^p depends on the product nia-uja, as shown in 
Fig. ISjb). This explains why the threshold m^/mi for the 
three-shell structure to occur in Fig.|3b) is roughly twice that 
in Fig. HJa). The (green) vertical dotted line indicates where 
the ^Li-^'^K mixture resides. In this case, a three-shell struc- 
ture appears in the upper panel while only a regular Sarma 
phase shows up in the lower panel of Fig. |5] 

It is worth pointing out that at large mass imbalance, the 
strong disparity between R^p and Rj^p makes population 
balance or imbalance less important. 

We end by noting that we have not included the species 
dependent, incoherent part of the fermion self energy, which 
may induce polarons in the mixed normal states, e.g., in the 
inner core of the three-shell structured phases. However, it 
is not important for the present study whether or not the 
minority species in these states form polarons. Following 
common practice 12. S IJs'-'^s'l, we have also neglected the 
particle-hole channel contributions [46], which can be roughly 
approximated by a shift in the pairing interaction strength 
147. 48^1 . These approximations are expected to modify the 
phase boundaries only quantitatively. In addition, we have 
not considered the FFLO states which, in an equal-mass Fermi 
gas, appears to be of less interest in 3D |13, i^fs^. Effects 
of mass imbalance on the FFLO phases will be investigated in 
a future work. 

In summary, we have studied the finite temperature phase 
diagrams for Fermi gases in a trap with both mass and pop- 
ulation imbalances, using a pairing fluctuation theory, with 
special attention paid to the ^Li-'^^K mixture. Unique to our 
theory are the wide spread pseudogap phenomena and the pre- 
diction of exotic phases, e.g., the phase-separated pseudogap 
phase, which can be tested by measuring the density and gap 
profiles in the trap. In particular, vortex measurements and rf 
spectroscopy |5ll] may be used to ascertain the superfluid and 
pseudogapped normal states. In order to compare with con- 
crete experiments, detailed parameters such as N^, uta, etc 
are needed. Our results can be tested experimentally when 
such experiments become available in the (near) future. 
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